Power-law rheology and mechano-sensing in a cytoskeleton 
model with forced protein unfolding 
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We describe a model of cytoskeletal mechanics based on the force- induced conformational change 
of protein cross-links in a stressed polymer network. Slow deformation of simulated networks con- 
taining cross- links that undergo repeated, serial domain unfolding leads to an unusual state — with 
many cross-links accumulating near the critical force for further unfolding. Thermal activation of 
these links gives rise to power-law rheology resembling the previously unexplained mechanical re- 
sponse of living cells. Moreover, we hypothesize that such protein cross-links function as biochemical 
mechano- sensors of cytoskeletal deformation. 
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The importance of mechanical cues for understanding 
cell behavior is increasingly recognized. Stem cell dif- 
ferentiation |l|, tissue morphogenesis 2] as well as cell 
growth Q and death Q are known to be affected by cell 
shape or the stiffness of the surrounding extra-cellular 
matrix. The molecular mechanisms by which such me- 
chanical cues produce biochemical responses, however, 
remain essentially unknown 5] . Similarly, the anomalous 
mechanical response of cells defies explanation despite 
the identification of the cytoskeleton's major structural 
constituents and considerable modeling effort 0, 0- 
Rheological measurements on living cells yield a fre- 
quency dependent shear modulus that scales as a weak 
power-law, \G*(lu)\ ~ cj^, with 0.1 < < 0.25, span- 
ning many decades of frequency 0, liol Hl| . Such a 
mechanical response is rather unusual 12], having only 
been observed previously in seemingly unrelated materi- 
als such as foams and pastes. While networks of purified 
filamentary biopolymers (e.g. F-actin) generally have 
a frequency-independent plateau elasticity [l3(, studies 
of networks contai ning the cross-link proteins filamin 
0, [3 or a-actinin |l6j have reproduced the weak power- 
law rheology as well as other aspects of the cell response, 
such as stress-induced stiffening. These results suggest 
that the power-law rheology of cells may be due to cross- 
link conformational change or unbinding. 

In this Letter, we propose a simple cytoskeletal archi- 
tecture that may explain both cell rheology and mechano- 
sensing, based on the force-induced serial unfolding of 
domains in protein cross-links. We show, using simu- 
lation, that shearing a simplified model network leads 
spontaneously to the accumulation of many cross-links 
at tensions on the cusp of unfolding. Thermally acti- 
vated unfolding of these cross-links readily reproduces 
cells' power-law frequency dependent shear modulus with 
physically plausible model parameters. Comparable 
models based on forced unbinding (e.g. of cross-links' 
actin binding domains from actin) rather than unfold- 
ing do not produce power-law rheology. Moreover, this 
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FIG. 1: Schematic representation of serial unfolding and the 
cytoskeleton. (a) As cross-links are extended, domains se- 
rially unfold with typical force vs. extension curves having 
abrupt transitions at a critical force, F c , corresponding states 
are labeled, (b) Our model of the cytoskeleton consisting of 
a network of generic semi-flexible polymers and partially un- 
folded, extensible cross-links. 



unusual, near-critical arrangement of cross-links is sug- 
gestive of a mechano-sensory function. We hypothesize 
that, by modulating the binding of signaling species, un- 
folding cross-link domains function as the fundamental 
biochemical transducers of cytoskeletal deformation. 

Many cross-linking proteins (e.g. filamin, a-actinin, 
spectrin, plakin, etc.) have a structure consisting of sim- 
ilar, repeated domains that can be serially unfolded by 
applied force [U 113 ? see Figure ^a). The unfolding 
time is typically described by the Bell model: |2£| 

t b (F) = r a eME B (l - y)\,E B = (1) 

where r a is a molecular attempt time, F c is a critical 
force, r Q is a characteristic bond length-scale, and fc^T is 
the thermal energy. The exponential form causes small 
changes in tension to yield rather large changes in un- 
folding time. As domains serially unfold, the molecule 
becomes progressively longer and longer. The entropic 
elasticity of the unfolded protein causes a spring-like re- 
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sponse between unfolding events, producing a 'sawtooth' 
force-extension profile, Figure [IJl 

We hypothesize that some cross-links unfold under 
physiological stresses, as sketched in Figure Gt, and that 
such unfolding is the source of microscopic stress re- 
laxation causing weak power-law rheology. Simple esti- 
mates |2l| suggest that physiologically plausible stresses 
of order ~100 Pa yield tensions sufficient to unfold some 
crosslink species. While the exponential sensitivity of 
the Bell model naturally gives rise to very broadly dis- 
tributed characteristic unfolding and relaxation times, a 
broad distribution alone is not sufficient. As we shall see, 
the power-law form of the rheology requires an exponen- 
tial distribution of cross-link tension. Just such force dis- 
tributions emerge spontaneously when simulated model 
networks are slowly sheared. 

For our simulations, we constructed two-dimensional 
networks having periodic connectivity in the left/right di- 
rection, whose nodes are initially offset from a triangular 
lattice by a Gaussian-distributed amount, Figure E{a,b). 
Force carrying mechanical links connect these nodes (15% 
of the nearest-neighbor bonds are left empty to introduce 
topological disorder). Networks larger than roughly 25 
by 35 nodes (~10 4 links) gave results that were essen- 
tially independent of system size and changes in network 
disorder. To model serial unfolding, each link has a linear 
'sawtooth' force-extension curve: F — k[{x— xq) mod x c ], 
where k is a spring constant, x is the instantaneous link 
extension, x c = F c /k is the critical extension for bond 
rupture and xq is the link length at the beginning of the 
simulation (typically 5-100 x c ). Our network does not 
contain separate rod-like and cross-link elements; each 
of our simulated links may be considered conceptually 
equivalent to a single inextensible rod of length xo con- 
nected in series with a cross-link having an infinite num- 
ber of unfoldable domains. Compressed links generate a 
compressive force (for x < x Ql F = —k[(x Q — x) mod x c ]), 
but this is not critical to our results. 

The network is sheared by translating the nodes 
clamped to the upper boundary in a series of small steps 
(A7 < 10 -3 ) at a constant strain rate, j(t) = jt. All 
nodes are first displaced according to an afflne shear 
deformation, and then relaxed to mechanical equilib- 
rium (zero total force on each node) by moving the non- 
boundary nodes using an 'over-damped' steepest descent 
algorithm. During relaxation, links move freely between 
branches of the force extension curve until an equilibrium 
configuration is reached. This is equivalent to assuming 
that domains unfold and refold instantaneously when the 
tension reaches F c and zero, respectively. 

To model Bell-type thermally-activated unfolding a 
Kinetic Monte Carlo (KMC) algorithm is used. The 
expected unfolding rate of each tensed (F > 0) link 
is computed from its tension and the Bell model (rate 
= 1/tb(F)). Consistent with the KMC algorithm, an 
exponentially distributed time step is generated which is 
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FIG. 2: Representations of our simplified model network 
and its force distribution without activated unfolding, (a) 
and (b) the network at strain 7 =0.3 and 1.5 respectively, 
(c) Probability distributions of scaled force, curves for strains 
7 = 0.3, 1.0, 1.5, and 2.0 (thin to thick) respectively. Line 
represents exponential behavior. 



inversely proportional to the total unfolding rate, and the 
link to unfold is selected in a rate weighted manner. To 
unfold the selected sub-critical link, its force extension 
curve is temporarily modified to be zero in the x c -wide 
extension interval it occupies. The network is then re- 
laxed, with the selected link typically moving into an 
adjacent, non-zero interval. 

It is illustrative to first consider the T = limit of 
the unfolding simulation (i.e. without Bell unfolding). 
In that case, the network configuration depends only on 
the strain, 7, and the network's geometrical parameters. 
For small strains, a roughly uniform distribution of link 
tension develops. For strains, 7 > 0.4, links accumu- 
late at tensions somewhat smaller than F c , Figure Etc). 
The distribution of link tensions approximates an in- 
creasing exponential function of force over a small range, 
P(F) ~ exp(F/F e ), the qualitative form required to pro- 
duce power-law rheology. As the strain is increased fur- 
ther, the range of forces showing exponential behavior 
broadens, but maintains the same slope: F e /F c w 0.4. 

In principle, the evolution of thermalized networks is 
more complicated, depending additionally on the Bell pa- 
rameters Eb and r a and the strain rate 7, (assumed to 
be slow, jr a < 1). In practice, however, the qualitative 
form of P(F) is quite similar to the athermal case, but 
with a somewhat smaller critical force, see Figure Et a )- 
Indeed, the P(F)'s from all of our thermalized unfolding 
simulations could be accurately scaled onto one another 
by renormalizing their forces by an effective critical force, 
F' c . Important for later discussion, the number of 'un- 
folded' domains in the network is a monotonic function 
of network strain, with little dependence on strain rate 
or hysteresis upon strain reversal, Figure Etb). 
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FIG. 3: Force distributions and scaling behavior with Bell- 
type activated unfolding, (a) Probability distributions of 
scaled force, determined for strain 7 = 1.5, Eb = 25 and 
r a = 50 nsec. Solid lines are from unfolding simulations with 
7 = 7, 1000 sec -1 and the athermal simulation (thin to thick) 
respectively. The dotted line is from the unbinding simula- 
tion at 7 = 0.5, Eb = 25, r a = 50 nsec, and 7 = 10 sec -1 (b) 
The number of unfolded domains as a function of strain. The 
solid is line is during strain application and the dotted line is 
during strain reversal, (c) Most probable unfolding force as a 
function of dimensionless strain rate. 



The effective critical force behavior has a simple expla- 
nation. Bell-type molecules subjected to a constant force 
loading rate have a well-defined most probable force for 
unfolding, which is a logarithmic function of the load- 
ing rate [l8|- The constant strain rate, 7, of our model 
slowly stretches links, leading to a constant force loading 
rate. We can identify the resulting most probable link 
unfolding force with F' c , which, as expected, scales loga- 
rithmically with 7, Figure Efc). Links with tensions even 
slightly smaller than F' c have an almost negligible Bell 
unfolding rate (since Eb ^> 1), leading P(F) to evolve 
much as an athermal model with 'all or nothing' unfold- 
ing at critical force F' c . 

The frequency-dependent mechanical response of our 
network can be estimated by superposing Maxwell modes 
[23. Formally, such a modulus describes the network's 
differential response to a small oscillatory applied stress 
superposed on the network's static stress. This is com- 
patible with the emerging view that cell measurements 
actually report such a differential shear modulus at a 
cell-generated prestress Q. As expected, the computed 
shear modulus, Figure Hftop), nas a roughly power-law 
form, |G*(o;)| ~ ujP, over a finite frequency range, with a 
non-universal exponent (3 ~ 2.bF c / 'F' c Eb- The modulus 
has a terminal mode at frequencies below the zero-force 
unfolding rate, r" 1 exp(— Eb), and approaches a plateau 
value, G at frequencies higher than uj ma x ~ 1/tb(F' c ). 
Surprisingly, the frequency range with power-law rhe- 
ology varies quite slowly with the strain rate: uJrnax ~ 
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FIG. 4: The computed frequency dependent mechanical re- 
sponse of our model, (top) Dimensionless differential stor- 
age (solid) and loss (dashed) moduli computed for Eb — 
25, 30, 35, 40 (thick to thin), (bottom) The experimental data 
of Fabry et al for untreated HASM cells, compared to our 
model (solid lines). At high frequencies, a G* ~ a; 0-75 , con- 
tribution was added (dashed line). Model rheology generated 
using Eb — 25, G Q = 5 kPa, r a = 70 nsec and 7 = 7 sec -1 . 
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In a separate study of forced unbinding, (rather than 
unfolding), force-extension curves with a single period 
were used (F = k(x — x Q ) for \F\ < F c , zero otherwise). 
Any link at zero force in mechanical equilibrium after a 
time/strain step corresponds to a link that has exceeded 
its critical force and unbound. We then reset its Xq pa _ 
rameter to its new equilibrium length to simulate rapid 
rebinding at zero force with the same network connec- 
tivity. Activated unbinding is handled as before, with 
the KMC-selected link having its force set to zero, fol- 
lowed by network relaxation and resetting of the x G 's of 
any unbound links. Unlike the unfolding simulations, the 
unbinding simulations never showed any accumulation of 
near critical links, FigureOfa), indicating such models do 
not produce power-law rheology. 

Before proceeding, some discussion of the model's 
strain rate, 7, is in order. We suppose that in cells 7 
is caused by molecular motor sliding or filament tread- 
milling at a roughly constant velocity. Rather than a 
uniform, pure shear deformation on the cellular scale, a 
spatially random deformation field with a typical, meso- 
scopic strain rate 7 would serve just as well. We fur- 
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ther suppose that other active processes continuously 're- 
model' the network with a turnover rate comparable to 
7, such that typical network segments are strained to 
7 ~ 1 in a dynamic steady state that replicates the non- 
steady behavior of our model at 7 ~ 1. To be consistent 
with 100-1000 second estimates for cytoskeletal turnover, 
7 should thus be in the range 10 -3 -10 -2 sec -1 . 

The computed rheology reproduces the cell response 
qualitatively, but falls short of quantitative replication. 
The experimental literature reports at least five fre- 
quency decades of power-law scaling with exponents in 
the range 0.10 < (3 < 0.25, while our model only 
yields five decades for (3 < 0.17 (or higher exponents 
over a narrower frequency range). This is due to the 
modest 'height' of our exponential force distribution, 
P(F' C )/P(F=Q) « 3, as in Figure EJa). Recent ather- 
mal simulations using a more realistic network structure 
yield a more pronounced exponential P(F) 23], suggest- 
ing our limited frequency range may be an artifact of 
our simplified network geometry. To estimate the model 
parameters corresponding to the physical case, we di- 
mensionalized our model, added a high-frequency con- 
tribution G* ^ cj 75 , and compared it to literature data 
. Our model reproduces the response of normal HASM 
cells, Figure bottom), with physically plausible param- 
eter values: Eb = 25, r a = 70 nsec and G Q = 5 kPa. The 
simulated strain rate, 7^7 sec -1 , is higher than our 
earlier estimate; smaller dimensionless strain rates were 
too computationally intensive to be accessible. Because 
of the slow variation of model rheology with 7, however, 
a more realistic value 7 = 10 -2 sec -1 , should yield a 
similar response with slightly smaller exponent, /3~0.15, 
and narrow the frequency range with power law rheology 
by (700)°- 18 = 3.2 x, or half a decade. 

The arrangement of structural molecules on the cusp 
of conformational change suggests an optimal configu- 
ration for a chemical mechano-sensor — suggesting that 
cells may maintain a metabolically costly dynamic cy- 
toskeleton as much for its sensory as its structural func- 
tions. Unlike other proposed sensor mechanisms 0, 13 
that transduce molecular stress, the earlier noted cor- 
respondence of unfolding and network strain indicates 
the sensing of deformation on the supramolecular scale. 
Many cross-linking species specifically bind a number of 
signaling proteins, including heat shock proteins, protein 
kinase C, Ral A J^IP2, PIP3, PI3-kinase, and MEKK1 
(for reviews [25l l2fij). Any of these proteins that specif- 
ically binds (or unbinds) cross-link domains upon forced 
unfolding would then transduce the shear or extensional 
strain of the network, which is presumably a prerequisite 
for shape or matrix compliance sensing. Furthermore, 
the myriad, multi-domain cross-link proteins and their 
isoforms localized to different parts of the cell suggest 
that each different cell sub-structure could have its own 
mechano-sensing capability. 

Our cross-link unfolding model reproduces the cell 



response and makes biochemically testable predictions: 
that some cross-link species should be partially unfolded 
under normal physiological conditions, and that their 
unfolding should increase with cell deformation. Such 
biochemical studies, combined with rheology simula- 
tions having more realistic network structure, hold the 
prospect of a cytoskeleton model grounded in polymer 
and single- molecule biophysics, one that can be inte- 
grated with mechano- sensory signaling pathways. 
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